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Abstract 

Examples of Morse functions with integrable gradient flows on some classical 
Riemannian manifolds are considered. In particular, we show that a generic height 
function on the natural embeddings of classical Lie groups and certain symmetric 
spaces is a perfect Morse function, i. e. has as many critical points as the homology 
requires, and the corresponding gradient flow can be described explicitly. This gives 
an explicit cell decomposition and geometric realization of the homology for such 
a manifold. As another application of the integrable Morse functions we give an 
(3JT)' elementary proof of Vassiljev's theorem on the flag join of Grassmannians. 

^ T + A + ■ 

^ , Introduction. 
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Morse theory based on a very natural idea till now remains one of the most attractive 
fields in the modern topology. We mention, for example, the papers by S. P. Novikov on a 
multivalued analogue of the Morse theory and those by A.T.Fomenko and his colloborators 
on the Morse theory for integrable Hamiltonian systems (see e. g. [|TJ). 

In the present paper we are not concerned with those generalizations but consider 
Morse theory in its classical aspects. It is well-known (see e. g. 0, 0) that it allows one to 
extract a topological information about a manifold M n from the knowledge of the critical 
points of some smooth function / : M n — > R. If one knows in addition the behaviour of 
the gradient flow of the function /, it provides a cell decomposition of the manifold M n . 

Note that for a given metric the gradient flow 

x = grad f(x) (0.1) 
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is a system of ordinary nonlinear differential equations, which in general does not admit an 
explicit integration. The critical points are defined by the system of nonlinear equations 
(in algebraic situation - by algebraic one) 

grad/(x) = 0, (0.2) 

which is also difficult to solve in general case. 

We will say that / is an integrable Morse function on a Riemannian manifold M 
if: first, / is a Morse function (or, more generally, Morse-Bott function) and, second, 
the equations ( p.l| ) and ( p.2| ) are "integrable". In the present paper we will mean the 
integrability in the naive sense as a possibility of explicit solution leaving apart more 
serious discussions of this notion for gradient flows. 

We should note that in the literature there exist already examples of gradient flows 
which could be naturally considered as integrable. For example, the Toda lattice can 
be interpreted as a gradient flow (see |]]). Once more remarkable example of an 
integrable gradient system was recently pointed out by S. P. Novikov in connection with 
the problem of "isoperiodical" deformations for finite-gap potentials and related problems 
of topological field theory || It has a feature that the corresponding metric is not 
positively definite. For discussions of the integrability of the algebraic equations ( p.2|) and 
examples of such systems we refer to 0, |[] . 

Such examples naturally appear in the theory of integrable systems. For instance, the 
stationary points of the mechanical system with a potential U(x) are nothing else but the 
critical points of U(x): 

grad U(x) = 0. 

It is natural to expect that for an integrable system the corresponding stationary problem 
also has an integrable nature. 

The starting point for us was another example arised from the theory of integrable 
systems with discrete time [U], |IT| . Such systems are defined by a "Lagrangian" 
which is a function £ on the product Q 2n = M n x M n . A sequence of points from M n 
X = {xk}, k G Z, is called a solution of the discrete system with Lagrangian £ if X is a 
critical point of the functional 

S(X) = Y,£(xk,x k+1 ). (0.3) 
kez 

In the papers || |10 , |Tl| it was shown that in the case when M n = SO(N) and 



£ = Tt(X JY l ), J = J 1 is a symmetric matrix, one has an integrable discrete version of the 
multidimensional top's dynamics. Such a form of the Lagrangian is uniquely determined 
by the following conditions: 

1) symmetry £(X,Y) = £(Y,X); 

2) left-invariance £(gX,gY) = £(X,Y), g G SO(N); 

3) bilinearity with respect to X and Y considered in the standard matrix realization. 
Note that because of the left-invariance £(X,Y) depends only on u = Y~ l X : 

£(X,Y) = £(Y- 1 X,Y~ 1 Y) = £(Y~ l X, I) = F(Y^X), where F is the corresponding 
function on the group, 

F(u) = Tr(Jcj). (0.4) 
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The critical points of the function F correspond to the stationary solutions of the discrete 
top's equations. 

It turns out that the function (p.4|) and its analogues on the unitary and symplectic 



groups [[ll]] can be considered as examples of integrable Morse functions on the correspond- 
ing group space with biinvariant (Cartan - Killing) metric. We discuss them in details in 
the first and second sections of the present paper. As an application we get a geometrical 



representation of the homology classes of these groups found by V. A. Vassiljev in fl2 |. 
For example, in the case of the unitary group we have the following isomorphism: 

n 

Hi{U(n)) * 0if^ 2 (G,(C n )), (0.5) 

k=0 

where Gk(C n ) is the Grassmannian manifold of complex fc-dimensional subspaces in C n . 

Let us define for an arbitrary Schubert cell a in Gk{C n ) the following cycle [a] in 
H 1r {U{n)). For each point in a, that is for a dimensional plane in C n , consider the set 
of all unitary operators preserving this plane and acting on its orthogonal complement as 



the reflection (— /) (in the paper [12] the trivial action on the complement is considered; 



this difference is not essential). The union of all such operators over all points from a 



forms a cycle [a]. Vassiljev's theorem |L2|, Theorem 1] states that the set of all cycles [a] 
corresponding to all Schubert cells of all Grassmannians Gk{C n ) gives a basis in the integer 
homology of the unitary group, which is reflected by equation ( |0.5| ). As V. A. Vassiljev 
told us this result was found also by M . Mahowald, so we will call equation (|0.5| ) Vassiljev- 
Mahowald decomposition. We will show that it can be obtained easily from the Morse 
theory for the function (p.4|) . 

We consider also the general height functions on symmetric embeddings of certain 
symmetric spaces like Grassmannians, Lagrangian Grassmannians, the spaces of complex 
and quaternian structures and show that they are perfect Morse functions with integrable 
gradient flows. For the Grassmannians the corresponding cell decomposition coincides 
with the classical Schubert decomposition. 

We should note that in a special case J = diag(l, 0, . . . , 0) function (|0.4|) was consid- 



ered and used for the investigation of the topology of classical Lie groups by L. S. Pontr- 



jagin in the paper [13 



It is interesting that the consideration of another integrable Morse function allows 
us to simplify essentially the proof of the following remarkable result from the same 
Vassiljev's paper fl2"fl . Let the Grassmannians Gi(V), G 2 (V),..., G n _i(V), where V is 
an n-dimensional vector space over k = R, C or quaternions H, be embedded into a 
Euclidean space of a large dimension N in a general position. The union 0„(k) of 
all the (n — 2)-dimensional simplexes with vertices (Vi, V2, ■ ■ ■ , V n -i), where the subspaces 
Vi G Gi{V) form a complete flag F, that is V\ C Vi C . . . C V n -i, endowed with the 
induced topology we will call the flag join of the Grassmannians. 



The theorem 5 from [T^ says that this space n (k) is homeomorphic to the sphere 
of certain dimension. In the simplest nontrival case 63 (R) can be identified with the 
quotient space of CP 2 by complex conjugation and, therefore, Vassiljev's result could 
be considered as higher-dimensional generalization of the Kuiper-Massey theorem stating 
that the last space is homeomorphic to S 4 . The proof given in [|T2J is based on some deep 
results of the modern topology. 
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In the third paragraph of the present paper we give an elementary proof of even a little 
more strong, smooth, version of this theorem. The idea is the following. Let's consider 
the set of the symmetric (Hermitian) matrices satisfying 

TrX =0 
TrX 2 = 1. 



This is a sphere of the dimension we need. Consider the function / = TrX 3 on it. Its 
critical points form in — 1) smooth submanifolds isomorphic to Gi(V), . . . , G?„_i(V). 
The corresponding gradient flow preserves the flag of the eigenspaces of a matrix X and 
determines in a natural way the decomposition of our sphere into "flag" simplexes a>, 
giving the equivalence with the flag join we need. It is interesting to note that the evolution 
of the eigenvalues of matrix X is described by the generalized Volterra chain (see. |T4| ). 

Our paper does not cover the whole variety of the applications of the integrable Morse 
functions. Looking at the history of the Morse theory from this point of view one can find 
a number of other remarkable examples. For instance, the well-known proof of the Bott 
periodicity JT|, |2| uses the functionals corresponding to integrable systems on the classical 
Lie groups. In the original Bott's proof it was the action functional on the loop space 
corresponding to the geodesic flow on the group which is obviously integrable, while in 
Fomenko's version [|15j the Dirichlet functional was used, whose extremals are known in the 
mathematical physics as the main chiral fields, the equations of which were integrated by 
the inverse scattering method only in 70-th (Zakharov and Mikhailov) jTjj, [I7|]. Of course, 
in both cases the integrability was used only to find explicitly some partial solutions, but 
anyway this also illustrates our approach to Morse theory. 

From this point of view it seems to be important that the equations of the famous 
self-dual Yang-Mills fields can be interpreted as the gradient flow of the Chern-Simons 
functional (see e.g. (T^|). It is well-known that with these equations are connected the 



recent remarkable achievements in topology (Donaldson, Witten, Seiberg). 

We note also an interesting paper by C. Tomei [|19| in which the Toda lattice is used 
for studying topology of manifold of isospectral Jacoby matrices. 



1. Height function on a symmetric space as a perfect Morse function. 

We start with consideration of the height functions on the classical Lie groups. Let G be 
one of the following compact Lie groups: 0(n), U(n), Spin). Here Spin) is the group of 
quaternian nxn matrices X such that X -1 = X . Sometimes it is called also the unitary 
symplectic group because it can be realized as the intersection of the complex symplectic 
group Sp(2n, C) with U(2n). 

Consider an arbitrary height function on G embedded in a natural way into the matrix 
space M{n, k), k = R, C or H, supplied with the standard scalar product 

(X,Y) = ReTr(X*F). 

Here for X £ M(rz, k) by X* we denote the matrix X . Such a function has the form: 

f A (X)=ReTr(AX), (1.1) 

for some matrix A e M{n, k). 
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Proposition 1.1 The critical points of the function f& are the matrices X G G that 
satisfy AX = (AX)*. 

The gradient flow for the function /a in the Cartan- Killing metric on the group G has 
the form 

X = A* - XAX. (1.2) 

Proof. Recall that the Cartan-Killing metric on the group G is the pullback of the 
Euclidean structure on the matrix space M(n, k) under the embedding. The right hand 
side of the equation flL2|) is the result of the projection of the vector A*, the gradient of 
the function in the ambient space M(n, k), onto the tangent space TxG: x{\[X~ l A*) — 

(X- l A*)*^j = A* - XAX, since X~ l = X* for X e G. The last expression vanishes iff 
AX = (AX)*.> 

Corollary 1.1 The gradient flow of the function Fa on the group G solves the problem 
of the polar decomposition for the matrix A. 

Indeed, let 

A = JQ, 

where J = J* is a positive definite matrix and Q e G, is the polar decomposition of the 
non-degenerate matrix A. Then for any X from G 

ReTr(JX) < ReTr(J) = ReTr(AQ*), 

so, X = Q* is the maximum point of the function Ja on G. It is not difficult to show (see 
also below) that the other critical points are not local maxima. Therefore, for almost all 
initial points X(0) = Xq e G for the equation (|1.2|) we will have X(t) — » Q* (t — > oo). 

Remark 1. It is not necessary to take as a starting point an element of the group G. For 
example, if we set X = then solving ( |1.2| ) we will always have X(t) — > Q* (t — > oo). 

To investigate the structure of the set of the critical points of the function J^, its 
gradient flow, and the corresponding cell decomposition of the group G it is sufficient to 
consider only diagonal positive matrices A = diag(a 1? . . . , a n ) since any nondegenerate 
matrix can be reduced to such a form by means of suitable left and right actions of G 
corresponding to some rotation of the ambient space M(n, k). 

Let's consider first the case when A — I, I is the identity matrix: fi(X) = ReTrX. 
It is easy to prove the following result (see for the details e.g. E( 



Proposition 1.2 1) The set of the critical points of the function fj on G consists of 
the involutions: X 2 = I and, therefore, can be naturally identified with the disjoint 
union of the Grassmannian manifolds G n)m (k), m = 0,1,..., n. Namely, to an 
arbitrary subspace V C k n is related the involutive operator Sy, the reflection with 
respect to V : 

Sv 



= I 

v 



S- 



V 



and all the involutions can be represented in such a way; 
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2) The separatrices of the gradient flow coming from the point S^eG fill the domain 
with closure isomorphic to 0(m), U(m), or Sp(m) correspondingly, m = dimV, 
that consists of the operators X such that X ^ ± = —I; 

3) the index of a critical point Sy, dim V = m, equals to 

{^m{m — 1) , k = R 
m 2 , k = C . 

m(2m + l) ,k = H 

Let's consider now the case when all the eigenvalues of the matrix A are different. 

Theorem 1.1 Let A = diaef £Zi , . . . , CL n ), < a\ < . . . < a n . Then 

1) $a is a Morse function on the group G; 

2) the critical points of the function f ' a are the matrices of the form 

X = diag(£i, . . . , e n ), where £j = ±f , i — 1, . . . , n; 

3) the index of such a critical point equals to 

it 

ind x (/ A ) = S> efeil ((dim R k)fc-l); (1.3) 

k=l 

4 ) the function / A has the minimal possible number of critical points among all Morse 
functions on G. 

Example. For the matrix 

A = diag(d — 1, 2c? — 1, . . . , nd — 1), d = diniR, k 

the value of the function f a at a critical point coincides with the signature of the Hessian 
at this point. This implies that the critical points of the same index lie on the same level 
and the greater index the greater level. Such functions are called sometimes Morse-Smale 
functions. 

Proof. We begin with the item 2). According to proposition 1.1 , if X is a critical point of 
f A then XAX = A = X*AX* = X^AX' 1 , therefore, A 2 = XA 2 X~ l and the matrices 
X and A 2 commute. It follows from the diagonal form of the matrix A that such matrix 
X should be diagonal too. 

Thus, X = diag(£i, . . . , e„), where £j = ±1. Let's compute the second derivative of 
the function f a at the point X . Consider the exponential map B i — ► X exp(_B), where 
B G TjG. The diagonal and off-diagonal elements of the matrix B can be taken as the 
coordinates in a small neighbourhood of the point X under that map. We have: 

f A (X) - f A (X ) = ReTr A(X - X ) = ReTr AX (exp B - I) = 

n 

= ReTr AX (B + ±B 2 ) + o(\\B\\ 2 ) = £ a^Ajbji + o(\\B\\ 2 ) = 
= - E (a i e i + a j e j )\b ij \ 2 + o(\\B\\ 2 ). 

n>j>i>l 
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Since dj > a,i > for j > i one has sgn(aj£j + dj£j) = sgn(cijEj). Thus \bij\ 2 , j > i, is 
contained in d 2 f with the negative sign if Bj = 1 and with the positive one if Sj = — 1. 
This implies the items 1) and 3). 

The item 4) follows from the Morse inequalities [[J valid for an arbitrary Morse function 

/: 

rrij > bj, (1.4) 

where rrij is the number of the critical points of index j of the function / and bj is the 
j-th Betti number. For the function J'a a simple calculation shows that the number rrij 
coincides with the corresponding Betti number (in the case G = 0(n) one should consider 
the homology groups over Z 2 ). The Betti numbers of the classical Lie groups have been 
calculated first by L. S. Pontrjagin in [13], one can find them also e. g. in |2l). Thus the 



Morse inequalities ( |1.4| ) for are in fact equalities and this completes the proof.> 

A Morse function with the minimal numbers of the critical points is called a perfect 
Morse function. We will give now a geometrical proof of this property for the height 
functions, which does not use the information about the homology groups of G in a more 
general situation. 

Let M be a symmetric space embedded in Euclidean Space E. Associate with each 
point x G M the reflection S x of the space E with respect to the normal subspace at the 
point x: 



= —I, dS x 

T X M 



We will call an embedding M — > E symmetric if all such reflections S x map M into 
itself. 

The classical groups G with the natural embeddings give the examples of such symmet- 
ric spaces: M = G, E = M(n, k), Sx{Y) = XY*X. Among other symmetric embedding^ 
we mention the following: 

1. Grassmann manifolds M = Gr njm (k) embedded in the space of symmetric (Her- 
mitian) matrices n x n as the set of the reflections with respect to the corresponding 
subspaces, Sx(Y) = XYX. 

2. Lagrangian Grassmannians, consisting of the Lagrangian subspaces in R 2n = C n : 
M = U(n)/0(n). It is embedded into the space of complex symmetric matrices as the 
subset of the unitary matrices, X = X T ,X~ l = X*. The reflection Sx has the form: 
S X (Y) =XYX. 

3. The symmetric space of complex structures in R 2n : M = 0(2n)/U(n). One can 
realise it as the set of the skewsymmetric orthogonal matrices in the Euclidean space of 
all skewsymmetric matrices, Sx{Y) = —XYX. 

4. The space of quaternian structures in C 2n : M = U(2n)/ Sp(n), realised as the 
set of the skewsymmetric unitary matrices in the space of all skewsymmetric matrices, 
S X (Y) = -XYX. 

5. Symmetric space M = Sp(n)/U(n) can be embedded into the space of skew- 
Hermitian quaternian matrices as the intersection with Sp(n), Sx{Y) = —XYX. 

The restriction to M of a linear function on E will be referred to as a height function. 
Among them the Morse functions form an open everywhere dense set. For example, 
the function f^ on the classical Lie group G is a Morse function iff the eigenvalues of 



1 All such embeddings are classified, see p8[ and Note added in proof below. 
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the positive symmetric (Hermitian) matrix J from the polar decomposition A = JQ are 
pairwise different (see the calculation above). 



Theorem 1.2 A Morse height function on a symmetric embedding of symmetric space 
M is a perfect Morse function. 

Proof. Notice first of all that by a small variation of such a function h which does 
not change the number and indeces of the critical points we can always reduce to the 
case when all the critical values of the function h are different and the gradient separatrix 
surfaces outcoming from the critical points intersect the incoming ones transversally. In 
particular, the separatrices outgoing from a critical point of index k never hit into critical 
the points of indeces > k and the set of such separatrices outgoing from a point of index 
k and incoming to a point of index {k — 1) is finite. 

Such a function h provides a cell decomposition of the manifold M in the following 
way. With each critical point of index k we associate a A;- dimensional cell which interior 
is the union of this critical point and all outgoing gradient separatrices. The incidence 
number of two cells corresponding to critical points of indeces k and (k — 1) is defined as 
the number of gradient separatrices outgoing from the first point and hitting the second 
one taken with some sign depending on the orientation of the cells, which is not essential 
for us. 

Let x be one of the critical points of h . Since dh(T X0 M) = the function h is 
invariant under the reflection S xo : S* Q h = h. Thus the set of the critical points of the 
function h _ _ is mapped to itself under such reflection. Since there is no more than one 



M 



critical point on the same level of h all the critical points are fixed under S. 



XQ ■ 



Consider now the cell complex C* defined by the function h . The groups Ck of this 
complex consist of the formal sums of the critical points of index k. Let points x and y 
have indeces k and (k — 1) correspondingly. As it was shown before, S x (y) = y. Since 
dS x = —I S x acts freely on the set of the gradient separatrices going from x to y and 

from S% — I it follows that the number of these trajectories is even. 

Thus all the incidence numbers of the complex are even and boundaries of all 
elements of the complex C* <E> Z 2 are 0. According to the Morse theory, the homology of 
this complex coincides with if*(M, Z2). The fact that the differential of the complex is 
trivial implies that Morse inequalities turn out to be equalities and the number of critical 
points of the function h is the minimal possible one among all Morse functions. > 

Remark 2. It turns out that for the height functions on the symmetric spaces the 
gradient separatrices never hit a critical point with "wrong" index, so we have a correct 
cell decomposition for all Morse height functions (see the next section). 

Remark 3. The height functions turn out to be perfect Morse functions not only for 
symmetric spaces but also for some other classical manifolds. For instance, this is the 
case for the Stiefel manifolds V n ,k naturally embedded into the space of n x k matrices 
and for the flag varieties realized as the set of symmetric (Hermitian) matrices with given 
spectrum. 
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2. Gradient flows and cell decompositions. 

The functions considered in the previous section have the following remarkable property: 
their gradient flows can be linearized by an appropriate change of variables and, therefore, 
can be integrated explicitly. 

Consider first the gradient flow of the function /a(X) = Tr AX, A = A* on the group 
G = 0(n), U(n), Sp(n): 

X = A-XAX. (2.1) 

Let 

y = (i - x)(i + xy 1 

be the Cayley transform of matrix X. It is known that the Cayley transform is involu- 
tive and provides a bijection between skew-symmetric (skew-hermitian) and orthogonal 
(unitary, symplectic) matrices not having —1 as an eigenvalue. 

Lemma 2.1 (|22j) The Cayley transform linearizes the flow (|2.1| ): 

Y = —(AY + YA). (2.2) 

Proof. 

Y = ((/ - X)(I + X)- 1 )' = -X(I + X)- 1 -(I- X)(I + X)- l X(I + X)- 1 = 
= -(/ + Y)X(I + X)- 1 =-\{I + Y)(A - XAX)(I + Y) = 
= + Y)(A- (I + Y)-\I -Y)A(I -Y)(I + Y)- 1 )(I + Y) = 

= + Y)A(I + Y)-(I- Y)A(I - Y)) = -(AY + YA).t> 



Proposition 2.1 The equation ( |2.1|) can be solved explicitly for arbitrary initial data 
X(0) = Xo- The solution has the form: 

X(t) = (sinh(At) + cosh(At)X )(cosh(At) + sinh(At)X ) _1 , (2.3) 

where sinh and cosh denote the standard hyperbolic functions of a matrix. 

Proof. For an initial matrix X not having —1 as an eigenvalue the equation ( |2.3|) 
is obtained by means of the Cayley transform from the solution of the equation fl2.2| ): 
Y(t) = exp(-At)Y exp(-At), with the initial point Y = (J - X )(I + X ) _1 . The set of 
the matrices not having —1 in the spectrum is everywhere dense in G, so the formula ( |2.3| ) 
gives the solution for all X it has a sense. 

The function thx = sinhx/ cosh a; takes values less than 1 for all real x, hence the 
matrix cosh(At) + sinh(At)X = cosh(At)(^I + th(^4i)XoJ is non-degenerate for all X 
from G.> 

Proposition 2.2 If the matrices A\ and A 2 are symmetric (Hermitian) and A\A 2 = 
A 2 Ai then the gradient flows of the functions f^ and f a 2 on the group G commute. 



9 



Proof. 



[grad f Al , grad f M \ \ x = (V grad f Al (grad f M ) - V grad fAa (grad f Al 
= V A ,-x Al x{ A 2 ~ XA 2 X) - Va 2 -xa 2 x(A 1 - XA X X) = 
= -(A x - XA 1 X)A 2 X - XA 2 (A l - XA-^X)+ 
+{A 2 - XA 2 X)A,X + XA t (A 2 - XA 2 X) = 
= (A 2 Ai - A ± A 2 )X + X{A X A 2 - A 2 A t ). 

The last expression vanishes for commuting matrices A\ and A 2 x> 

Corollary 2.1 Grassmann manifolds G n}Tn (k) embedded into G are invariant under the 
gradient flow (|2.1|) . The restricted flow coincides with the gradient flow of the function 



f A with respect to the induced metric on G n m (k). 

(k) X ......V / 

Proof. Since the gradient flows of functions // and f A commute the last flow should 
preserve the set of the first one's stable points. Note that G ntTn (k) is a connected compo- 



nent of the critical point set of the function // (see. proposition |1.2|) . The restricted flow 
coincides with the gradient flow of the function f A 
metric on G n m (k) because of the following general result. > 



with respect to the induced 



Proposition 2.3 Suppose that a submanifold N of a Riemannian manifold M is invari- 
ant under the gradient flow of some function f on M. Then the restricted flow on N 
is the gradient one with respect to the restricted function f 
metric on N. 



and induced Riemannian 

N 



Proof. It is easy to see that in general case the gradient flow of the restricted function 
/ coinsides with the orthogonal projection of grad / in the induced metric to the tagent 
subspace of iV.> 

It is easy to check that all the symmetric spaces listed in the previous section are invari- 
ant under the gradient flows of the height functions on the corresponding Lie groups. In- 
deed, for the Lagrangian Grassmannians LG n = U (n)/0(n) this follows immediately from 
the flow equation ( |2.1| ). For the spaces of comlex and quaternian structures 0(2n)/U(n) 
and U(2n)/Sp(n) one can see it from the equations ( |1.2|) : 



X = A* - XAX = -(A + XAX) = -X T 
since A T = -A X T = -X. Similarly, for spaces Sp(n)/U(n) A* = -A, X* = -X, 

X = A* - XAX = -(A + XAX) = -X*. 



Theorem 2.1 The gradient flows of the height functions on the classical Lie groups and 
the embedded symmetric spaces described above can be integrated explicitly. 
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Proof. Using an appropriate shift by an orthogonal (unitary, symplectic) matrix one can 
always reduce the equations of the gradient flow to the form ( |2.1|) and then apply the 
explicit formula fl2.3|) .> 

For the Grassmannians one can give another explicit formula for the gradient flow. 
Let us determine an element V £ Gr„ iTO (k) bynxm matrix Z with the columns forming a 
basis in the subspace V. Then the operator of the reflection with respect to the subspace 
V has the form: 

X = 2Z(Z*Z)- 1 Z* - I. 

Lemma 2.2 If Z(t) is changing according to the equation 

Z = AZ, 

then X(t) satisfies (|2.1| ). 
Proof. Indeed, 

X = (2Z(Z*Z)- 1 Z* - /)' = 

= 2Z(Z*Z)- 1 Z* + 2Z(Z*Z)- 1 Z* - 2Z(Z*Z)- 1 (Z*Z + z*z)(z*zy 1 z* = 
= 2AZ(Z*Z)~ 1 Z* + 2Z(Z*Z)- 1 Z*A - 4Z(Z*Z)~ 1 Z*AZ(Z*Z)- 1 Z* = 
= A- XAX.> 



Proposition 2.4 The evolution of a subspace V £ G ntm (k) under the gradient flow of 
the function f a on the Grassmannian Gr njm (k) is given by 

V(t) = exp(At)V. 

Let's consider now the cell decompositions determined by our Morse functions. 
Let A = diag(a 1; . . . , a n ), < a x < a 2 < . . . < a n . As a simple corollary of the the 
previous description of the gradient flow one has 

Proposition 2.5 The cell decomposition defined by the function }a on the Grassmannian 
G n|7n (k) coincides with the classical Schubert decomposition. 

The following statement describes the corresponding cell decomposition of the group 



G, which turns out to be a small modification of the one considered by V. Vassiljev [12 . 



Theorem 2.2 The cells of the classical Lie group G determined by the Morse function 
f a correspond bijectively to the Schubert cells of Grassmannians G nj?Ti (k) ; m = 0, 1, . . . , n. 
The interior of the cell corresponding to a Schubert cell a consists of the operators such 
that the eigenspace corresponding to the eigenvalue —1 has the orthogonal complement 
belonging to a. 
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Proof. As follows from the proposition LI all the critical points of the function lie on 
the Grassmannians embedded into G. This gives the bijectivity. It is easy to check that 
the dimensions of the cells coincide with the indeces of the corresponding critical points. 
Thus only the invariance of the cells under the gradient flow has to be proven. Let X(t) 
be a gradient trajectory. Denote by V_(i) the eigenspace of the operator X(t) with the 
eigenvalue —1. We state that 



which implies 



V4t)=exp(-At)(V(0)), 
(V.(t)) 1 = e X p(At)(y(0)) ± . 

(cosh(At) - sinh(Ai) 



V-(0) 



VI (0) 



Indeed, 

(cosh(At) + sinh(At)X(0) 

This means that 

(cosh(At) + sinh(At)X(0))~ 1 (exp(-At)V_(0)) = VL(0). 
Now taking into account (|2.3|) we get: 
X(t) 

= ( sinh(At) + cosh(At)X(0)) ( cosh(At) + sinh(At)X(0) N 



exp(-At) 



(V-(0))' 



(sinh(At) - cosh(v4t)) exp(At) 



cx P (-At)y_(o) 

exp(-At)V-(0)' 



cxp(-At)V-(0) 

Thus exp(-At)(V(0)) C V(t). In a similar way C exp(-At)(V(0)). This com- 
pletes the proof.> 

We see that the cell decomposition is common for a big family of the Morse functions 
J'a- This is a particular case of the following general result. 

Theorem 2.3 Any two functions with commuting gradient flows from a connected fam- 
ily of Morse functions on a compact Riemannian manifold M determine the same cell 
decomposition of M. 

Proof. From the commutativity it follows that the gradient flow of the first function 
preserves the stationary set of the second gradient flow. But since they are Morse functions 
this is a discrete set and, therefore, both functions have the same critical points. For given 
point x the cell containing this point corresponds to the critical point which is a limit at 
the plus infinity of the gradient trajectory starting from x (see above). It is clear that for 
functions close enough we have the same limit and, therefore, the same cell decomposition. 
Now the theorem follows from the connectness of the family.> 

This allows one to know the cell decomposition for such a family if one knows it for a 
particular function from it. For example, in the family we can choose 

A = diag(d — 1, 2c? — 1, . . . , nd — 1), d — diniR, k 

when the corresponding f'A is a Morse-Smale function (see section 1). 

As a corollary of the previous results we have the Vassiljev-Mahowald decomposition. 
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Theorem 2.4 ( ||12|| ) For the homology of the classical Lie groups one has the following 
isomorphisms: 

H l (0(n),Z 2 ) S ®H. Kk-i)(G n . k (R),Z 2 ) 
Hi(U(n),Z) = ©iZ i _ fc2 (G rM ,(C),Z) 

fc=0 
n 

Hi(Sp(n),Z) - © fl i _ fc(2fc+1) (G n , fc (H),Z). 

fc=0 

In a similar way the following result can be proved. 

Theorem 2.5 For t/ie homology of Lagrangian Grassmannians, spaces of complex and 
quaternian structures and spaces Spin) /U{n) one has the following decompositions: 

n 

H i (U(n)/0(n),Z 2 ) = H fc(fc+1) (G w . fe (R), Z 2 ) 

n 

Hi(0(2n)/U(n),Z) = ^_ fc(fc _ 1) (G„, fc (C), Z) 

fc=0 

H i {U{2n)/Sp{n),Z) S © ^_ 2fc2 (G n , fc (H), Z) 

fc=0 

n 

Hi(3p(n)/U(n),Z) = H i _ fc(fc+1) (G n , fc (C), Z). 

fc=0 

For Lagrangian Grassmannians the cell decomposition determined by a suitable height 
function turns out to coincide with one sometimes called Schubert decomposition.^ It 
could be described as follows (see [p3fl): for given Schubert cell a in G n( jt(R) we define 
the cell [a] C LG n as the set of all Lagrangian subspaces L in C n with a standard 
symplectic structure such that the projection V of L to R n C C n belong to a. Indeed, 
let U = X + iY be a unitary matrix, which columns represent a basis in L. One has 
(X + iY)(X T - iY T ) = I and, therefore, XX T + YY T = I and YX T - XY T = 0. This 
implies that the symmetric unitary matrix UU T acts as — I on the orthogonal complement 
of V because UU T = (XX T - YY T ) + i(YX T + XY T ) and X T v = for all veV ± . 

Cell decomposition of the space Sp(n)/U(n) can be constructed in a way analogous 
to Latrangian Grassmanians. Geometrically the space Spin) /U in) can be realized as 
the set of complex n-dimensional subspaces in H n , on which the restriction of bilinear 
skewsymmetric function Q(u,v) = (uj,v) is trivial, where (■,■) is the hermitian scalar 
product in C 2n . To any such a subspace L a matrix U = X + jY from Spin), where 
X, Y e M(n, C), corresponds columns of which make up a basis in L. Matrix UiU* is 
skew-Hermitian (as a quaternian one) and the operator —iUilf* acts as —I on V- 1 -, where 
V C C n is the subspace spanned over C by the columns of matrix X. 

For the space of complex structures CS n the gradient cell decomposition implying 
isomorphism mentioned above is constructed in the following way. Let's fix some complex 
structure Jo £ CSn, that is an orthogonal operator in R 2n , satisfying the relation Jq = — /. 
Then for a Schubert cell a C G nj fe(C) we define [a] as the union of all complex structures 

2 We would like to quote D. B. Fuks |2j| : "The manifold U(n) /0(n) has a convenient cell decomposition 
which also could be referred to as Schubert one although it was invented by Arnold." 
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J preserving V and coinsiding with J on V 1 - for all V G a. In the quaternian case we 
have a similar construction. 

We would like to mention that one can prove other decomposition formulas using a 
suitable Morse-Bott height functions with critical submanifolds. For instanse, using the 
gradient flow of the height function f& with diagonal A having only two different eigenval- 
ues we can prove the following decomposition formula for the homology of Grassannians. 

Theorem 2.6 The following isomorphisms hold for arbitrary n = n% + 

Hi(G nt k(R), Z2) — © flj_( m _n, 1 )fe a (G ni( fc 1 (R) x G n2) jfc 2 (R), Z 2 ) 

ki+k2=k 

Hi(G ntk (C),Z) = © Hi_ 2 

(ni — fci)fe (^ni,fei 

(C) X G n2i A; 2 (C),Z) 

^(G? n , fc (H),Z) = © Hi — 4(ni — k\)k2 (^nijfei 

(H)xG n2jfe2 (H),Z). 



A flow similar to ours was used by V. M. Buchstaber |23], Lemma 3.7] in theory of 
characteristic classes. 

3. Flag join of Grassmannians and an integrable gradient flow on the sphere 

In this section we will give an elementary proof of the following geometrical result of 
V. Vassiljev. Let's consider a general embedding of the Grassmannians Cr n>1 (k), G nj2 (k), 
. . . , G n .„_i(k) into a Euclidean space of a large dimension. The union of the simplexes 
with the vertices V\ € G n> i(k), V2 G (j n> 2(k),. . . , V n -\ G G n>n _i(k) forming a complete 
flag: Vi C V2 C . . . C Kj-i, with the induced topology we will call the flag join of 
Grassmannians and denote following to Vassiljev by n (k). 



Theorem 3.1 (V. Vassiljev ||12| . Theorem 5]) The flag join of Grassmannians is iso- 



morphic to the sphere of the appropriate dimension: 

e„(R) ^ 5^+- 2 

n (C) ~ S n(n-l)+n-2 
9 n (H) ~ S 2n(n-l)+n-2^ 

To prove this let's consider the function / = | Tr X 3 on the set of symmetric (Hermi- 
tian) n x n matrices V satisfying the conditions: 



TrX = 
TrX 2 = 1 



2 _ 1 (3.1) 



This is the sphere S N of dimension 

N = n-2 + (dim R k) n ^~ 1 \ 
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Proposition 3.1 The set of the critical points of the function f on the sphere S N is 
a disjoint union of the Grassmann manifolds G n ,i(k), G n ^(\s), ■ ■ ■ , G n ,n-i(k) smoothly 
embedded into S N in such a way that to a subspace V G G njm (k) corresponds the matrix 
X = X(V) such that 

X 
X 



I 

V V nm 

(3.2) 

I. 



y± ' V n ( n ~ m ) 

Proof. Functions TrX, TrX 2 , TrX 3 in the space of symmetric (Hermitian) matrices 
have gradients: 

grad(TrX) = I 

grad(TrX 2 ) = IX 

grad(TrX 3 ) = 3X 2 . 

The critical points of the function / are those X for which matrices /, X, X 2 are 
dependent. This holds iff X has no more than two different eigenvalues. Since TrX = 0, 
TrX 2 = 1 such a matrix X has exactly two different eigenvalues \i\ and /i2, /ii < < ji2- 
Let m and (n — m) be their multiplicity. It is easy to find from ( |3.1| ) that 



n — m I m 



nm y n(n — m) 

An operator with such a spectrum is uniquely determined by its eigenspace corresponding 
to n\.> 

Lemma 3.1 The gradient flow of the function f has the form: 

X = X 2 - f(X)X --I. (3.3) 

n 

Proof. The right hand side of the equation ( |3.3|) is the orthogonal projection of the 
vector X 2 to the tangent space of our S with the normal vectors X and /: 

X 2 - <^4a' - !£4l = X 2 - ^lx - ™I/ = X 2 - (Tt X>)X - -I. 
(X,X) (/,/) TrX 2 TrI K ' n 

It is well-known that it gives the gradient vector of the restricted function with respect 
to the induced metrio 

Proposition 3.2 The gradient flow of the function f preserves the eigenspaces and the 
anharmonic ratios of the eigenvalues of the matrix X . 

Proof. The first statement follows from the fact that the left hand side of the equal- 
ity (|3.3 ) is a polynomial of X with scalar coefficients. The dynamics of the eigenvalues 
is determined by a vector field of the form (a\ 2 + bX + c)Jv, where (a, 6, c) ^ (0,0,0), 
which is the same for all of them. But it is well-known that such a vector field generates 
a one-parametric group of the projective transformations of the line.> 
Let 

Ai(X)<A 2 (X)<...<A n (X) (3.4) 



15 



be the eigenvalues of the matrix X, 



V^X), V 2 (X),..., V n (X) 
be the corresponding eigenspaces. Let's define the eigenflag of X as 

F(X) = {U Q CU 1 C...C U n }, (3.5) 

where 

k 

U (X) = {0}, U k (X) = ®Vi(X), l<k<n. 

i=l 

If some eigenvalues are equal to each other the corresponding eigenflag will be incomplete. 



Proposition [T2] gives the set of integrals of the flow ( |3.3| ), consisting of the eigenflag 
of X and the anharmonic ratios of the eigenvalues. This set is complete in the sense that 
if we fix the values of the integrals we will have a trajectory. 

Indeed, a symmetric (Hermitian) matrix is completely determined by its spectrum 
and the eigenflag. Given the anharmonic ratios the spectrum is defined up to a projective 
transformation. The group of projective transforms of the line is 3-dimensional, therefore, 
the conditions 



n 



E Ai = 

i= n l (3.6) 

EAf = 1, 

i=l 



equivalent to ( |3.1| ) determine a 1-dimensional family. 

For two flags F\ and F 2 we will write F\ < F 2 if all the subspaces of F\ are contained 
also in F 2 . 

Lemma 3.2 For a given flag F = {U = {0} C Ui C U 2 C . . . C U k+ i = k™} the 
matrices X satisfying fl3.1| ) and such that F(X) < F form a (k — 1)- dimensional simplex 
with the vertices corresponding to the subspaces Ui, ...,£/& according to (|3.2|). The ratios 



_ A m (X)-A t (X) 

l ~ \ n {x)-\,{x) ' {6J) 

are the barycentric coordinates on this simplex. 

Proof. For simplicity we give the proof only for the case of complete flag F. In 
this case we can assume without loss of generality that the subspaces Ux,...,U n are 
<ei>, <ei, e 2 >, . . . , <e\, . . . , e n >, where ei, . . . , e n is the standard basis in k n . Then the 
condition F(X) < F is equivalent to the statement that X is a diagonal real matrix with 
the increasing elements at the diagonal: X n < X 22 < . . . < X nn . 

So, there is one-to-one correspondence between the set of matrices X e S N satisfying 
F(X) < F and the set of sequences Ai < . . . < X n , satisfying (|3.1|). We should only check 
that the functions (|3.7| ) can be considered as the barycentric coordinates on this set. 

Indeed, by construction we have 

71—1 

J2 a i = 

i=l 

a,; > for all i — 1 , . . . , n — 1. 
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The sequence {Aj} is denned by {c^} from Q3.7|) uniquely up to an affine transform. The 
normalization conditions ( |3.6| ) together with the increasing property ( |3.4| ) determine it 
uniquely. > 

Summarizing, we have the following 

Theorem 3.2 The gradient flow of the function f = \trX z on the sphere in the space 
of symmetric matrices X : trX = 0, trX 2 = 1, is integrable and gives a decomposition 
of this sphere into a union of simplexes. The vertices of such a simplex belong to the 
Grassmannians embedded into the sphere as the critical set of f and form a flag. 

Concerning the integrability we would like to mention that the gradient flow of the 
function / in the barycentric variables a% after a suitable change of time has the form of 



the generalised Volterra chain (see [|14 

d 

—ai = di 
dr 



i— 1 i+n 

ah~ 53 a i I ' ( 3 - 8 ) 



\k=i— n l=i+l 



where we set with i < and i > n to be zero. 
Indeed, from (13. 31) we have: 



1 n 

\i = A 2 - J 3 A 4 , where h = Y] A 3 . 

n fc=i 



Hence 



^ f K \ _ K+i — Ai (A n — Ai)(A^ + i — Xj) _ I X i+1 — \j \ n — Ai \ _ 

1 \ X n — X\ ) X n — X\ (A„, — Ai) 2 * yAj + i — Aj X n — Ai J 

= ( X 2 l+1 - hX l+l - (A 2 - IjX) _ A 2 ,-/ 3 A W -(A 2 -/ 3 A 1 ) \ = 
1 \ Aj + i — Aj A n — Ai J 

(i-l n-1 \ 

k=l l=i+l J 

Let's change now the time variable on the gradient flow trajectory: dr = (A n — Xi)dt. 
The equation will take the form: 



d_ 

ar 



t— 1 n— 1 



Ys a k- a i I ' ( 3 - 9 ) 

fc=l l=i+l 



which coincides with ( |3.8| ) if one set ao = ct-i = a-2 — ■ ■ ■ — a n — a n+i = a-n+2 = . . . = 0. 

In this case the generalized Volterra chain can be easily explicitly integrated. To do 
this let us change variables (cf. 

, Aj + i — Ai 

hi = a { 



k=l 



A n — Ai 



We will have: 



d , , Xi + \ + Ai — X n — Ai 

-i-h = b, = bi{bi - 1), 

dr X n - Ai 
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hence 



where the constants Cj are determined by the initial data: q = 1 — 6j(0) . The last 
formula is in a good agreement with our previous considerations. 

As a corollary we have the following strengthened version of Vassiljev's theorem. 

Theorem 3.3 Flag join of the Grassmannians n (k) admits the smooth structure of the 
standard sphere for which all the embeddings G„, TO ^ ©n(k) are smooth. 

Acknowledgements. We would like to thank V. M. Buchstaber for pointing out the 
references [[H], [24]. One of us (A. V.) is grateful to J. Moser and P. Santini for useful 
discussions and to Forschungsinstitut fur Mathematik (ETH, Zurich) for the hospitality 
during the summer term of 1995. 



Note added in proof. Recently L. Nicolaescu attracted our attention to the papers 
25| , |26| , f27f , where the properties of the height Morse functions on a class of so-called R- 



symmetric spaces are investigated. These spaces turn out to be precisely the ones for which 
embeddings with the property we demanded in section 1 do exist |2l| . It is interesting 
that the same embeddings are remarkable from another geometrical point of view: they 



are minimal in the sense of total curvature (see p9| , |30|). The whole list of R-symmetric 



spaces besides the examples we mentioned in section 1 consists of M = Sp(n)/U(n) 
and exceptional spaces M = E§/ Spin(10) x T, F±j Spin(9) or Cayley projective plane, 
E-j/Eq x T and Eq/F^ (see f29fl). It would be interesting to investigate the integrability 
of the gradient flows for the height functions on the corresponding embeddings. We are 
very grateful to L. Nicolaescu and E.Leuzinger for pointing out these remarkable papers 
to us. 
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